Faraday Rotation and Circular Dichroism Spectra of Gold and Silver Nanoparticle 
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We study the magneto-optical response of noble metal nanoparticle clusters. We consider the 
interaction between the light-induced dipoles of particles. In the presence of a magnetic field, the 
simplest achiral cluster, a dimer, exhibits circular dichroism (CD). The CD of a dimer depends on 
the directions of the magnetic field and the light wave vector. The CD of a populous cluster weakly 
depends on the magnetic field. Upon scattering from the cluster, an incident linearly polarized light 
with polarization azimuth ip, becomes elliptically polarized. The polarization azimuth rotation and 
ellipticity angle variation are sinusoidal functions of 2<p. The anisotropy and the chirality of the 
cluster control the amplitude and offset of these sinusoidal functions. The Faraday rotation and 
Faraday ellipticity are also sinusoidal functions of 2tp. Near the surface plasmon frequency, Faraday 
rotation and Faraday ellipticity increase. 

PACS numbers: 78.20.Ls,78.67.Bf,73.20.Mf,78.20.Fm 
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I. INTRODUCTION 

The magneto-optical (MO) response of nanocompos- 
ites has attracted a lot of interest. In view of po- 
tential applications in information storage, magnetic 
field sensors, switches, modulators, optical isolators, cir- 
culators, and deflectors^, composites containing 7 — 
Fe2C>3, CoFe204, Fe304, Co, and Fe nanoparticles 
are studied^—. Recent advances in plasmonics have 
led to the introduction of magnetoplasmonic nanos- 
tructures composed of ferromagnets and noble met- 
als. The large MO response of magnetoplasmonic com- 
posites can be understood heuristically: The metal- 
lic constituent sustains localized or propagating sur- 
face plasmons. This gives rise to enhanced electro- 
magnetic field intensity inside the ferromagnetic con- 
stituent which exhibits the MO activity. Magnetoplas- 
monic crystals 1 ^, metal films with arrays of holes filled 
with a magneto-optically active materia l 11 ' 12 , Fe/Cu bi- 
layer nlms±2, Ag/Co/Ag and Au/Fe/Au system s 14 ' 15 , 
Au-coated 'y-Fe20s nanoparticle s 16 ' 17 , and dumbbell- like 
Ag-CoFe2 04 nanoparticle pair s 18 ' 19 are thoroughly stud- 
ied. 

Recently, Sepulveda et al. 20 have experimentally 
demonstrated that pure noble metal nanoparticles ex- 
hibit measurable MO response at low magnetic fields 
(about 10 4 G). This effect was attributed to the Lorentz 
force induced by the collective motion of the conduction 
electrons in the nanoparticle when the localized surface 
plasmon resonance is excited. In this experiment, the 
particles were not in close proximity, thus interparticle 
coupling effect was weak. 

In this paper, we study interparticle coupling effect on 
the MO activity of gold and silver nanoparticle aggre- 
gates. We rely on the celebrated dipole approximation. 
In this approximation, a particle with diameter much less 
than the light wavelength is represented as a dipole. The 
local electric field acting on any dipole is a superposi- 
tion of the incident light field and secondary fields pro- 



duced by other dipoles. We assume that the aggregate 
interacts with the external static magnetic field B and 
linearly polarized light specified by the polarization az- 
imuth if. Using the azimuth ip s and ellipticity angle fi 
to characterize the vibration ellipse of the field scattered 
by the aggregate, we find that the polarization azimuth 
rotation A = ip s — ip and ellipticity angle variation fj, can 
be expressed as 



Afa) = A st » + BA Pa », 



(1) 



The term A atr (/z atr ) which does not depend on the mag- 
netic field, describes the structural polarization azimuth 
rotation (structural ellipticity angle variation), a mani- 
festation of the anisotropy and chirality of the aggregate. 
The term -BA Far (_B/i Far ) which is proportional to the 
magnetic field, describes the Faraday rotation (Faraday 
ellipticity). We find that near the plasmon resonance of 
the noble metal, both A Far and /x Far increase. We find 
that 



A„ tr [>) pa Aisin(2^-2^ ) + A 2 , 
Mrtrfa) ~ Mi sin(2<£ - 2<p' ) + M2, 



(2) 



where Aj and /i^ depend on the light frequency, host 
medium dielectric constant, nanoparticle radius, etc.. 
A Far and [i Far are also sinusoidal function of 2<p. For 
the simplest anisotropic and achiral cluster, a dimer, we 
find that A2 = and /j,2 = 0. For a set of nanoparticles 
positioned on a helix, A2 7^ and fi2 7^ 0. For a random 
gas of particles, which is anisotropic and chiral to some 
extent, A2 7^ and /12 7^ 0. The significance of these re- 
sults may be seen better by considering optical response 
of a bulk medium. For a bulk medium, the azimuth rota- 
tion is also a sinusoidal function of 2tp. The amplitude Ai 
and the offset A2 are largely controlled by the anisotropy 
and chirality of the bulk medium, respectively^ 2 -. 

Optical rotation and circular dichroism (CD) spectra 
of bulk medium are not independent, but are closely 
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connected by Kramers-Kronig relations^ 3 -. For metal 
nanoparticle assemblies, this interrelationship is not yet 
thoroughly studied. Drachev et alM> observed large local 
CD in fractal aggregates of silver nanoparticles by means 
of photon scanning tunneling microscopy. The value of 
CD averaged over the whole sample was close to zero. 
Yannopapas 2 ^ noted that interacting nanospheres occu- 
pying the sites of a rectangular lattice, show a consid- 
erable CD around the surface plasmon frequency. This 
planar nonchiral structure is indeed extrinsically chiral 26 . 
Fan and Govoro v 27 ' 28 showed that the helical arrange- 
ments of gold nanoparticles exhibit CD in the visible 
wavelength range. It is noteworthy that chiral molecules 
and biomolecules show strong CD in the UV range. Cir- 
cular dichroism of metallic nanoparticle assemblies orig- 
inates in the electromagnetic dipole-dipole interaction 
between particles. Recently Kuzyk et a/.— used DNA 
origami to create left- and right-handed nanohelices of 
diameter 34 nm and helical pitch 57 nm, by nine gold 
nanoparticles of radius 5 nm. 

Inspired by the work of Sepiilveda et al.—, and Fan 
and Govoro v 27 ' 28 , we study CD of interacting metallic 
nanoparticles subject to an external static magnetic field. 
We find that the simplest achiral cluster, a dimer, ex- 
hibits CD in the presence of the magnetic field. Near 
the surface plasmon frequency, CD increase considerably. 
Quite remarkably, CD depends on the directions of the 
magnetic field B and the light wave vector k with re- 
spect to the dimer axis n. In particular, CD vanishes 
if (i) k || n and B _L n, (ii) k _L n and B \\ n, and 
(iii) k _L n, B _L n, and k _L B. This reminds us the 
intriguing magnetochiral dichroism 3 ^— of a solution of 
chiral molecules: The absorption of a light beam parallel 
to the magnetic field and a beam antiparallel to the field 
are slightly different. Indeed, on the basis of symmetry 
arguments, it has been shown^ that the dielectric tensor 
of a chiral medium has a term proportional to k ■ B. 

Our article is organized as follows. In Sec. [TT] we in- 
troduce the model. The circular dichroism and Faraday 
rotation spectra of a dimer, a set of nanoparticles posi- 
tioned on a helix, and a random gas of nanoparticles are 
presented in Sees. IIII1 IIV1 and E] respectively. A sum- 
mary of our results and discussions are given in Sec. I VII 



II. MODEL 

A. Polarizability tensor of the nanoparticle 

We adopt the free electron Drude model to describe the 
magneto-optical response of bulk metals. For an electron 
subject to the static magnetic field B = B(b x ,b y , b z ) and 
the oscillating electric field E(t) — Re[E(ui) exp(— iuit)] of 
the light, the equation of motion reads 




Here x, 7, m e , and — e < denote the displacement vec- 
tor, relaxation constant, effective mass, and charge of the 
electron, respectively. CGS units are used throughout 
this paper. 

The polarization of the medium is P{t) = —en e x(t) 
where n e is the number of conduction electrons per 
unit volume. Solving Eq. ([3]) with the ansatz 
x(t) = Rc[x(u) exp(— iut)} yields P(to) = —en e x(uj) — 
x{u)E(uj). For weak magnetic fields (about 10 4 G) the 
cyclotron frequency 

^ = ^- (4) 

m e c 

is much smaller than the light frequency u. Thus we 
calculate the electric susceptibility tensor and the 

permittivity tensor e(uS) = I + Attx(uj) to first order in 
uJc/u. Indeed 

(s xx (u>) \Gui c b z -\Guj c b y \ 
-\Guj c b z e yy (uj) \Guj c b x , (5) 
iGuj c by -iGu c b x e zz {w) ) 

where e xx = e yy = e zz = 1 - w 2 /(w 2 + vyu), G = 
w 2 /(w 2 + i7w) 2 , and uj p = y' '4irn e e 2 /m e is the plasma 
frequency. 

Now we consider a metal nanoparticle of radius a much 
smaller than the light wavelength. The polarizability ten- 
sor of the spherical nanoparticle is^ 

a(u) =a 3 {e-E m I)(e + 2s m I)-\ (6) 

where e m is the dielectric constant of the isotropic host 
medium and J is the identity matrix. To first order in 

(a xx (uj) -Tu c b z Tu c b v \ 
Tu c b z a yy (u!) -Fuj c b x , (7) 
-Fio c by F(jj c b x a zz (uj) J 

where a xx = a yy = a zz = a 3 (e xx -e m )/ (e xx + 2s m ) and 



+ i7) 2 (e xx + 2e m ) 2 ' 

The off-diagonal terms of a which are proportional to 
the magnetic field B, might seem small at first glance. 
However, the factor J- becomes considerably large at fre- 
quencies where Re(e xx + 2e m ) = ±Im(e xx ). These fre- 
quencies are near to the surface plasmon frequency where 
Re(e xx + 2e m ) = 0. 

The simple Drude model helps us to understand the 
origin of plasmon-induced magneto-optical activity of no- 
bel metal clusters^?. To make a better link to exper- 
iments, we use the Drude critical points model— for 
the dielectric function of bulk Ag and Au. This model 
well reproduces the experimental data of Johnson and 
Christy 38 . We also modify the dielectric function of 
bulk metal to take into account the finite size effects in 
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small metallic spheres. We rely on the following radius- 
dependent dielectric function^ 



_bulk 



W + ITbulk^ U + 17W 



(9) 



where 7 = Tbulk + Av F /a is the radius-dependent relax- 
ation, vp is the Fermi velocity, and the parameter A is 
of order 1. To calculate the size corrected a xx and J 7 , 
we assume that A = 3/4^S, although there is a slight 
disagreement among authors about A^— . For silver, 

and 



lo.„ = 1.386 x 10 



16 „-l 



7buik = 4.584 x 10 irf s 



13 „-l 



16 „-l 



v F = 1.39 x 10 s cm/s. For gold, uj p = 1.32 x 10 lb s 
7buik = 1-08 x 10 14 s" 1 , and v F = 1.4 x 10 8 cm/s. 



B. Coupled-dipole equations 

To describe interaction of an incident electromagnetic 
wave with a set of N identical particles, we use the dipole 
approximation: Each metal nanoparticle of radius a be- 
haves as a point dipole with polarizability tensor a. The 
local field acting on any dipole is a superposition of the 
incident field and secondary fields produced by other 
dipoles. Thus the coupled-dipole equations (CDE) for 
the induced dipoles are 

N 

d l {u J )^a{u){E{r l ,u)+Y J W ^- r j)- d A^))^) 

i=i 

Here r< and Re[di(u>) exp(— iuit)} denote the position and 
light-induced dipole of the ith particle, respectively. The 
dyadic Green's function W determines the field produced 
by the oscillating dipole. 



W a ,/s(r«)=[(l-iA!|r <J -|)( 



2 S a t 



£ m\fij\ 



e ifc l r «l,(ll) 



where r%j = — rj, k = wyf^/c, and the Greek sub- 
scripts denote the Cartesian components. 



III. DIMER 

A. Circular dichroism 

It is instructive to consider a dimer, as the sim- 
plest achiral cluster. We assume that the particles re- 
side at r\ = (0,0,0) and r 2 = (L, 0, 0). The vector 
n = (r 2 — T\)/L specifies the axis of the dimer. The 
dimer interacts with the electromagnetic field E(r,ui) — 
Eoeexp(ik ■ r). The polarization vector e and the wave 
vector k = fc(sin#cos(/3,sin#sin<y9,cos#) are perpendicu- 
lar. 



The light-induced dipole moments obey the equations 
d 1 =a(E(r 1 ,u)+W-d 2 ^ and d 2 = a\E(r 2 , u)+W -dA . 



Here 



W 



'77 + £ N 
T) 
77/ 



1 \h p ikL 



L e m L 
3 3i/c s e ifeL 



£(L,k) = i-k z + - 1 -—) 



L 2 



(12) 



Introducing d m = a 1 (di — d 2 ), d p = a 1 (di + d 2 ), 
E m = E(n,,u>) - E(r 2 ,u), and E p = E(n, w) + 
E(r 2 ,ui), the coupled-dipole equations can be written as 

(J + Wot)d m = E m , 

(I-Wa)d p =E P , (13) 

which are easily solvable. To first order in lu c /lu, 

_(l^a xx )E mx -(rt+£)a xy E my -(rt+£)a xz E mz 



(1 + r/a xx )(l + {r) + £)a xx ) 



1 / YjOLxy^mx ^j^yz^mz 

my — T~ I ~ ' 7 \ 7\ r^mij Z ' 

l-tyctxx V1 + (77 + Qol xx ±+r}a xx 



1 



rja xz Er, 



r\a. yz E. 



yz-^iny 



-En 



l+V^xx VI + (r/ + S)a xx l+rja xx 
(l-r]a xx )E px +(r]+£)a X yE P y+(rr+£)a xz E pz 



dpX (1 -??Q!xx)(l - {V + £)°<*x) 

. _ 1 ( mxyEpx jp ya y zE pz 
\-r\a xx \ l-{rj + i)a xx l-i]a xx J 



d„ 



1 



7]a xz E p 



rja yz E p 

l-Wzz V 1 - (r/ + i)a xx l-rja x . 



■+E pz ).(U) 



The absorption cross section of the dimer can be ex- 
pressed in terms of the auxiliary variables d m and d p . 
Indeed 



Air to 



hn{cT l *dl ■ d l + oT l *d* 2 ■ d 2 ) 



-^0 \J C 

= W^HF- • d m + ad P ■ d *p)- ( 15 ) 

V m 

We pay attention to the absorption cross sections a± 
corresponding to the polarization vectors 

e± = —= (cos 8 cos tp =F i sin ip, cos 9 sin <p ± i cos ip, — sin 9) . 
V2 

e + and e describe the right and left circularly polar- 
ized light, respectively. A straightforward but lengthy 
calculation shows that the dimer is circular dichroic. 

<7+-(7_ = lo c [k x b x Kx (77, £)+(k y by+k z b z )A 2 (r], £)J k 2 x L 2 

+ u c ( K k x b x k 1 {-j 1 ,-i)+{kyby+k z b z )k 2 {-ri,-£,)) 
x(4 + k 2 x L 2 ), (17) 
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FIG. 1: (<7+— CT-)/(7ra 2 ) of an Ag dimer as a function of the 
wavelength 2tvc/oj for various directions of B and k. Here 
a =10 nm, L = 3a, e m = 1.7 and B = 8000 G. 



where 



47T 



47T 



Rc(^)+Im( i 27?J " )Im(tt M ) 



-Im(- 



+(rj + £)a xx 



S m \l+mxx\ 2 



-Re(- 



ll-KrH^axxl 1 



)Re(J") 
)Im(a :c2 



(18) 



As expected, the achiral dimer exhibits no circular 
dichroism in the absence of the external magnetic field 
(u c = 0). Quite remarkably, c+ — <r_ depends on the 
directions of the magnetic field B = B(b x ,b y ,b z ) and the 
wave vector k = (k x ,k y ,k z ) with respect to the dimer axis 
n, see Fig. [T] The circular dichroism certainly vanishes 
if (i) k || n and Bin, (ii) k _L n and B \\ n, and (iii) 
fc 1 n, B 1 n, and fc _L £?. 

Note that Ai oc J and A2 oc T . Near the surface 
plasmon frequency, T and hence <r+ — er_ increase con- 
siderably. Moreover, one may expect enhanced circular 
dichroism where Ke(l±rja xx ) = or Ke(l±(r]+^)a xx ) = 0. 
We found that these resonances hardly occur for inter- 
particle distance L > 3a. Note that the dipole approxi- 
mation is not reliable in the region L < 3a. 

Figure []] exemplifies er+ — <r_ of an Ag dimer for vari- 
ous directions of B and k. Here a = 10 nm, L = 3a, and 
B = 8000 G. For fe||(0,0, 1) and B||(1,0,0) the dimer ex- 
hibits no circular dichroism, but upon light propagation 
along fc|| (0.5,0,0.866) the dimer responds quite differ- 
ently. Considering the case fc||(0, 0, 1) and -B||(0, 0, 1), we 
find that the peaks of c+ — cr_ may increase and shift as 
the direction of magnetic field varies. 

Figure [5] demonstrates the circular dichroism <r + — <r_ 
in units of the geometric cross section 7ra 2 , as a function 
of the wavelength 2itc/uj. Here a = 10 nm, k\\ z and 
B = 8000 Gz. An Ag dimer exhibits a considerable neg- 
ative (positive) dichroism at wavelength 2ttc/uj — 383 nm 
(391.5 nm) when L = 3a and e m — 1.7, see Fig. HJa). 
These peaks increase, and slightly shift to lower wave- 
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FIG. 2: (a+-a-)/{na 2 ) as a function of the wavelength 2nc/u) 
for (a) Ag dimer in a medium with e m = 1.7, (b) Au dimer in 
a medium with e m = 1.7, and (c) Ag dimer with interparticle 
distance L = 3a. Here a —10 nm, k\\z and £? = 8000 Gz. 



lengths as L increases. The spectrum does not change 
significantly upon increasing the interparticle distance 
beyond L = 15a, since the dipolar interaction becomes 
quite weak. In the same host medium, the dichroism 
of an Au dimer reaches its maximum at wavelength 
2ttc/ui = 540 nm, see Fig. &h). The dichroism of Ag 
dimer is about 20 times stronger than that of Au dimer. 
This is expected, as l/lm(e xx ) and hence T of Ag is 
greater than that of Au. Figures [2lc) shows the impact 
of host medium on the circular dichroism spectrum of 
Ag dimer. Both peaks of the spectrum grow and shift to 
higher wavelengths as e m increases. 



B. Faraday rotation 

Now we assume that the dimer interacts with a linearly 
polarized light with the polarization vector 



(cos ip, sin tp, 0). 



(19) 
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The wave vector k and the external magnetic field B are 
along the unit vector z. We measure the electric field E s 
radiated by the dimer at (0,0, h), where h ^> L > a. In 
the far (radiation) zone the scattered field takes on the 
form 



(d\ — zz ■ d\ 



2Akh 



k 2 e 



(ad p — zz ■ adp) 



k 2 e ikh 



(d 2 -zz ■ d 2 ) 



k 2 e 



e m h 



Evnh 



(20) 



The electromagnetic wave E s = (E sx , E sy , 0) is ellip- 
tically polarized. The vibration ellipse can be character- 
ized by its azimuth ip s , the angle between the semimajor 
axis and the unit vector x, and its ellipticity tan/x, the 
ratio of the length of its semiminor axis to that of its 
semimajor axis. Operationally defined in terms of mea- 
surable intensities, the Stokes parameter a 43 ! 44 



I 

Q 
u 

V 



\E SX \ 2 + \E sy \ 2 , 

\E gX \ 2 - \E S y\ , 

2Re(E* sx E sy ), 
2Im(E* sx E sy ), 



(21) 



determine completely the state of polarization. Indeed 



1 U 
<p s = - arctan—, 

1 . V 
H = — arcsm — . 



(22) 



Note that the azimuth and the ellipticity angle of the 
linearly polarized incident wave are tp and 0, respectively. 
Hereafter we focus on the polarization azimuth rotation 
<ys s — <ys and the ellipticity angle variation fj,. 

Powered by the analytical expressions (fT4|) and (|20| . 
we find that to first order in uj c /uj 



1 Jisin2<ys 

= —arctan — w 

2 v cos2^+J2 ; 

w c J\Jz sin 2 2ip+ ( J3COS 2ip 
+ 1> (Ji sin2<p) 2 + (cc 



~ r ,-J4)(cos2</5+ J 2 ) 
(Ji sin2^) 2 + (cos2v?+ J 2 ) 2 



1 . . — J5 sin 2ip 
A* = « arcsml — ) 

2 V 2 cos2(p+r 

; (Je cos 2ip+Jr)(J 2 cos2</?+l)— J4J5 sin 2 2(p 

(J 2 cos 2ip+l)y/(J 2 cos 2(p+l) 2 - ( J 5 sin 2^) 2 

Co 



V ~2 



(23) 



where the parameters J1-J7 are reported in the appendix. 
As already mentioned in the Introduction, the polariza- 
tion azimuth rotation is composed of two terms which 
are different in nature. The first term of A = ip s — ip 
which does not depend on the external magnetic field, 
is the structural polarization azimuth rotation A otr (</?). 
Note that a dimer is an anisotropic and achiral object. 
Switching off the dipolar interaction between the spheres 
(?7,£ — > 0), this structural azimuth rotation expectedly 
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FIG. 3: A 5tr ((/p) of (a) Ag dimer, (b) Au dimer, and (c) Ag 
dimer, as a function of tp for various wavelengths 2nc/oj. Here 
a = 10 nm, L = 3a, e m = 1-7, and B — 0. 



vanishes. The second term of <p s —(p which is proportional 
to the external magnetic field, is the Faraday rotation 
B A Far (ip) . Note that J3 oc T and J4 oc J 7 , thus the plas- 
mon resonance boosts the Faraday rotation. Similarly, 
the first term of \i describes the structural ellipticity an- 
gle variation. J4, Jq and J7 are proportional to J- . The 
plasmon excitation thus enhances the second term of /i 
which is the Faraday ellipticity 

We have numerically studied J1-J7. We find that in a 
large domain of parameter space 

Jl~l, J 2 <1, J 3 W C <1, J 4 W C «1, J 5 <1, J 6 Wc<l, J7^ C <1, 

(24) 

which allows us to rewrite Eq. (|23|) as 



<^ s - ^3 sin 2^9+ ^ ( {Ji - J2J3) cos 2y -I- ./..). 



— sin 2tp + ^y(— J?J 2 cos 2ip + J 7 ). (25) 



/' 



In other words, provided that condition (|24|) is satisfied, 
A 5tr , i?A Far , jti atr and Bfi Fai are sinusoidal function of 
2ip. These functions can be represented in the form of 
Eq. @: For a dimer A 2 = and \x 2 = 0. 

Now we use Eq. (|2Uf to access the structural polariza- 
tion azimuth rotation A . of a silver dimer as a function 
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FIG. 4: BA Far (y>) of (a) Ag dimer, (b) Au dimer, and (c) 
Ag dimer, as a function of ip for various wavelengths 2ttc/lj. 
Here a =10 nm, L = 3a, e ro = 1.7, and B = 8000 G. 
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FIG. 5: (a) ^i str (<£>) of Ag dimer, (b) Bfj, FaT (ip) of Ag dimer, 
and (c) B/i Far {p>) of Au dimer, as a function of <p for various 
wavelengths 2nc/u). Here a = 10 nm, L = 3a, e m = 1.7, and 
B = 8000 G. 



of <p and 2ttc/ui. We assume that a = 10 nm, L = 3a, 
e m = 1.7, and B = 0. Figure E[a) shows that A atr is a 
sinusoidal function of 2y>. At wavelength 420 nm, A str 
has a considerable amplitude 7.84°. Upon increasing the 
wavelength to 500 nm, the amplitude decreases to 3.62°. 
Figure O^b) shows a similar plot for a gold dimer. Here 
A atr is less than 4.1°. Figure [3fc) shows that in the nar- 
row window 380 ^ 2irc/u> ^ 395 nm, A str of a silver 
dimer is a nonsinusoidal function. This is because the 
condition (|24p is violated. Quite remarkably, in this win- 
dow A atr can be as large as 22.6°. FiguresSJa) and[3Jh) 
show the Faraday rotation of a silver and a gold dimer 
subject to the magnetic field B = 8000 G, respectively. 
BA FtLi ((p) is a sinusoidal function of 2ip. Apparently the 
maximum of Faraday rotation < 0.01° is much less than 
the maximum of structural azimuth rotation 7.84°. Fig- 
ure [4jc) shows that in the narrow window 380 — 400 nm, 
BA Fax (<p) of a silver dimer is a nonsinusoidal function. 

The structural ellipticity angle variation /i str behaves 
similar to A str . The sinusoidal /i atr of a silver (gold) dimer 
gains its maximum 0.1° (3.3°) at wavelength 480 nm 
(510 nm). In the window 380 — 395 nm, fj, str of a sil- 
ver dimer is a nonsinusoidal function, see Fig. [SJa). For 
the magnetic field B = 8000 G and wavelengths not in 



the window 330 — 440 nm (540 — 556 nm), S/i Far of a 
silver (gold) dimer is a sinusoidal function of 2ip. In this 
regime, the maximum of Faraday ellipticity angle varia- 
tion 0.004° (0.0046°) of a silver (gold) dimer is much less 
than the maximum of structural ellipticity angle varia- 
tion 0.1° (3.3°). The nonsinusoidal behavior of B/i Far for 
silver and gold dimer are shown in Figs. [5jb) and[5jc), 
respectively. 



IV. NANOPARTICLES POSITIONED ON A 
HELIX 

We aim to further study correlation between the chiral- 
ity and the magneto-optical response of a set of nanopar- 
ticles. Here we consider a set of N identical nanoparticles 
of radius a positioned on a helix. We assume that the ith 
particle resides at 

Xi = -Rhclix COs(i#hclix — 0hclix), 
Di = ^hclix sin(i6>holix — #hclix), 
Zi = (i - I)phelix0helix/(27T). (26) 
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FIG. 6: Magneto-optical response of 15 silver nanoparticles 
positioned on a helix, a = 10 nm and e m = 1.7. (a) (cr+ — 
<7_)/(7ra 2 ) as a function of 2-kc/uj, for various directions of 
k — fc(sin 9 cos 95, sin #sin <p, cos 6 1 ). Here _B = 0. (b) A Btr as a 
function of <p, for various wavelengths 2-kc/uj. (c) A2, J12, A4, 
and /I4 as a function of 2-kc/uj. Here _B = 8000 G. 
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FIG. 7: Magneto-optical response of 1000 silver nanoparti- 
cles distributed in an ellipsoid. a = 10 nm and e m = 1.7. (a) 
(<r_|_— cr_)/(7ro 2 ) as a function of 2-kc/uj. (b) A atr as a function 
of if, for various wavelengths 2nc/u. (c) A2, ~p2, A4, and 7Z4 
as a function of 2-kc/uj. Here fc || z and i? = 8000 Gz. 



As an example, we choose N = 15, a= 10 nm, -Rhciix = 
30 nm, phciix = 30 nm, and Rhciix = 7r/3. 

We assume that the helical structure interacts with 
the electromagnetic field E(r,u) — Eoeexp(ik ■ r). We 
compute the absorption cross sections o± corresponding 
to the polarization vectors e± introduced in Eq. (IT6l) . 
The dielectric constant of the matrix e m = 1.7 and ex- 
ternal magnetic field B = 0. Figure Uta) delineates 
<7_l_ — cr_ of a helical assembly of silver nanoparticles as 
a function of wavelength 2ttc/uj, for various directions of 
k = k(s'm6 cosip, smOs'mip, cos 8). We can see that the 
CD spectrum has both positive and negative bands. The 
peaks of the CD spectrum change and shift as the direc- 
tion of light wave vector k varies. For light propagation 
along the helix axis, i.e. k \\ z, <J+— cr_ gains its con- 
siderable extrema 6.057ra 2 and — 5.637ra 2 at wavelengths 
389 nm and 396 nm, respectively. We also find that^ the 
CD spectrum does not change significantly as the static 
magnetic field increases to B = 8000 G. We observe that 
the CD spectrum is quite susceptible to neglect of finite 
size effects in small metallic spheres, see Eq. ©. This 



point is overlooked in Ref. 

We now assume that k \\ z, B \\ z, and the helical 
structure interacts with a linearly polarized light with 
the polarization vector introduced in Eq. (jTTJ)) . Direct 
numerical solution of the coupled-dipole equations (fT0|) 
allows us to obtain the azimuth ip s and ellipticity angle fi 
which characterize the vibration ellipse of the scattered 
field. Figure EJb) shows the structural polarization az- 
imuth rotation A Btr (iy9) of a helical structure, as a func- 
tion of ip. We observe that A atr is a sinusoidal function of 
2ip. We also find thai— n str , A Far , and /i Far are sinusoidal 
functions of 2ip. Remarkably, for a set of nanoparticles 
positioned on a helix A2 ^ and \xi ^ 0, see Eq. ©. 

For an arbitrary set of N nanoparticles, we can not 
provide exact analytical expressions for A(<p) and n{<p), 
or numerically explore the whole parameter space. In 
other words, the (narrow) windows of parameters where 
A.(ip) and fi(ip) are nonsinusoidal function of 2<p, are not 
known. However, encouraged by the success of Eq. ([2) in 
describing the optical behavior of a bulk medium^ 2 " and a 
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dimer, we propose to measure the well denned quantities 

a; = r 

JO T 

- r r ^ 

Jo ^ 

si = r b^m-, 

Jo n 
TU = [ B^XV) — , (27) 

to study correlation between the geometry and the 
magneto-optical response of a set of nanoparticles. Fig- 
ure [6jc) shows A2, /IJ, A4 and pj of our helical structure 
as a function of the wavelength 2itc/uj. A 2 gains its min- 
imum —0.023° and its maximum 0.011° at wavelengths 
384 nm and 397 nm, respectively. JX2 reaches its maxi- 
mum 0.026° at wavelength 388 nm. Note that near the 
same wavelengths, the circular dichroism spectrum is ex- 
tremum. Similar to A2, A4 has a deep minimum —0.023° 
at 384 nm. However A2 and A4 are distinct. Notably, 
A2, JtJ) A4 and 7Z4 have the same order of magnitude. 



V. RANDOM GAS OF NANOPARTICLES 

We now consider a random arrangement of TV identical 
nanoparticles of radius a. We rely on Na 3 /Rg as the 
volume fraction of the cluster. The radius of gyration 
R g which characterizes the size of the system, is defined 

as R 2 g = jj J2iLi I r i. - r cm | 2 - Here r cm is the position 
vector of the center of mass. 

Intuitively, a random gas of nanoparticles is anisotropic 
and chiral to some extent. One can quantify chiral- 
ity of a set of points through the Hausdorff chirality 
measure^. To characterize the anisotropy of the cluster, 
one can calculate the moment of inertia tensor — 
J2iLi da^Ti ■ ri — (r i ) Q (r i ) ( g and its three real eigenval- 
ues I\ ^ I2 ^ I3. Then an appropriate function of the 
principal moments of inertia, e.g. + h + 13) serves 

as a measure of anisotropy. 

We study a set of 1000 silver nanoparticles of radius 
a = 10 nm, which are distributed in an ellipsoid with 
semi-axes lengths 10 /mi, 30 /mi, and 58 /mi. Deliber- 
ately, we did not allowed the particles to be nearer than 
3a. The volume fraction of the sample is about 0.05. We 
further assume that fe|| z, B || z and e m = 1.7. Our cluster 
exhibits circular dichroism, see Fig. [TJa). The spectrum 
has three maxima and three minima, whose position and 
magnitude are almost independent of B. We observe 
that^ A Btr , /i str , A Far , and /i Far are sinusoidal functions 
of 2ip, see Fig. EJb). Remarkably, for a random gas of 
nanoparticles all A2, JL2, A4 and /Z4 are nonzero, and 
have the same order of magnitude, see Fig. Etc). As ex- 
pected, the wavelengths at which er+ — <t_, A 2 , /tj, A 4 
and J14 are extremum, are quite close. 



VI. DISCUSSION 

On propagation through a bulk medium, a linearly po- 
larized electromagnetic wave becomes elliptically polar- 
ized. The polarization azimuth rotation and ellipticity 
angle variation are due to optical anisotropy and optical 
activity of the bulk medium. The constitutive equation 
AirP a — (e a p — S a p)Ep + r Q ^ 7 V 7 -Eg gives the polariza- 
tion of the bulk medium P in terms of the electric field 
E. Here e a p are components of the dielectric tensor de- 
scribing anisotropy and T a ^ are components of the non- 
locality tensor describing optical activity^ 2 -. In isotropic 
optically active media e a p = n 2 S a p and T a p 1 = Te a p 7 , 
where n is the complex refractive index and e Q ^ 7 is the 
Levi-Civita symbol. If a left (right) circularly polarized 
light is incident on the isotropic optically active medium, 
it will propagate as a left (right) circularly polarized light. 
The eigenpolarizations of an anisotropic medium with no 
optical activity (r Q( g 7 = 0) are two linearly polarized 
waves. In general, one can measure the polarization az- 
imuth rotation A and ellipticity angle variation /1 as a 
function of the azimuth (p of an initially polarized wave, 
to separate the optical anisotropy and the optical activity 
of bulk medium 2 ^. 

A nanoparticle assembly is anisotropic and chiral to 
some extent. We showed that, to a good approxima- 
tion, the structural polarization azimuth rotation and the 
structural ellipticity angle variation of a cluster are sinu- 
soidal functions of 2<p, see Eq. Studying the dimer, 
the helical arrangement and the random gas of nanopar- 
ticles, we inferred that the amplitude Ai and the offset 
A2 inform us about the anisotropy and chirality of the 
cluster, respectively. Further studies are needed to quan- 
tify anisotropy and chirality of the nanoparticle assembly 
and quantify their link to Ai and A2. Ai and A2 of a 
random gas of 1000 silver nanoparticles are of the order 
of 1° and 0.01°, respectively, see Figs.Etb) andEKc). Ex- 
perimental measurement of these angles is not difficult. 

In the presence of an external magnetic field 10000 G, 
bulk gold exhibits very weak MO response. Indeed at 
wavelength 700 nm, e xy = 10 _4 + il0 -3 for gold, whereas 
£xy = 0.8 + i0. 16 for cobalt 2 ^. However, the off-diagonal 
terms of the polarizability tensor of an isolated metal 
nanoparticle becomes considerably large near the surface 
plasmon frequency^ ., see Eqs. ([7]) and ©. Hui and 
Stroud studied Faraday rotation of a dilute suspension 
of particles in a host. Relying on the Maxwell- Garnett 
approximation, they also found that the Faraday effect 
becomes large near the surface plasmon frequency^. The 
effective medium or "mean- field" approximations do not 
take into account the interaction between particles. Thus 
we relied on coupled-dipole equations. Considering in- 
teracting nanoparticles, we found that Faraday rotation 
and Faraday ellipticity are also sinusoidal functions of 
2ip. Remarkably, Faraday rotation and Faraday elliptic- 
ity of both helical arrangement and random gas of silver 
nanoparticles are of the order 0.01°, see Figs. [He) and 

ETc). 
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Recently Kuzyk et al£2- used DNA origami to create 
left- and right-handed nanohelices of diameter 34 nm and 
helical pitch 57 nm, by nine gold nanoparticles of radius 
5 nm. They even achieved spectral tuning of CD by 
metal composition: Depositing several nanometers of sil- 
ver on gold nanoparticles, a blueshift of the CD peak was 
observed. In their pioneering work, Fan and Govorov 27 
assumed that helical arrangement of nanoparticles are 
in a solution and have random orientations. Thus they 
performed averaging over the various directions of light 
wave vector k to compute the CD spectrum. We en- 
visage an ensemble of aligned helical structures, thus we 
emphasized the dependence of the CD spectrum on the 
direction of k, see Fig. Eta). First, we remind that var- 
ious techniques are developed to align nanorods^, e.g. 
gold nanorods are incorporated within poly(vinyl alco- 
hol) thin films and subsequently aligned by heating and 
stretching the film 49 . Nanohelices maybe aligned by the 
same procedures. Second, we believe that the light pres- 
sure can be invoked to generate and align helical arrange- 
ment of nanoparticles. Nahal et al, 50 reported generation 
of spontaneous periodic structures in AgCl-Ag films. The 
AgCl layer acts as a waveguide. The incident light scat- 
tered by the Ag nanoparticles excites waveguide modes, 
which interfere with the incident beam. Silver nanopar- 
ticles then immigrate to the minima of the interference 
pattern. Electron microphotographs of samples irradi- 
ated by a circularly polarized light, suggests the forma- 
tion of helical arrangement of nanoparticles. This is fur- 
ther supported by the circular dichroism of the sample, 
and recent AFM and SEM images^. 

We anticipated that the simplest cluster, a dimer, 
has interesting MO properties. Two dimensional ar- 
rangements of dimers can be produced by electron beam 
lithography 52 . Molecular bridges can also be employed 
to form Au and Ag dimers, trimers, and tetramer a 53 ' 54 . 

Our work can be extended in many directions. MO 
response of trimers, tetramers, two-dimensional arrange- 
ments of particles, fractal clusters, and the influence of 



particle size dispersion^— , are of immediate interest. 
We envisage the optical control of the MO response of a 
solution of molecularly bridged dimers. Some molecular 
bridges, for example azobenzene derivatives, reversibly 
transform between trans and cis conformations upon 
light radiation. Thus the dimer length and its MO prop- 
erties can be optically controlled. 



Acknowledgments 

We would like to thank A. Nahal from University of 
Tehran, for fruitful discussions. 



Appendix A: The parameters J1-J7 



The parameters J1-J7 introduced in Eq. (|23[) can be 
written as 



Jo = \a xx -i]a 2 xx \ 2 + \a xx - (?7 + £) a L| 2 j 
2 

Ji = —Rc[{a xx -r]al x ){a* xx - (jl'+fKJ], 
Jo 



1 
2 

Jo 
2 



(kx-^Ll 2 - -(v + C) a Ll 2 )> 

Re[(a xx - V a 2 xx )T* + (a xx - (rf + fK^J] 



Jo 
2 



Re[(a xx - va 2 xx )T* - (a* x - W H*W XX )F] , 



J 2 
J3 
Ji 
J 5 

J& = -r^[( a xx-ml x )^*+( a *xx- (v* + rKiW. 

•Jo 

_2 

J 7 =—lm[(a xx -r 1 al x )J J *-(a xx -(if+C)a x 2 x )J : }. (Al) 
Jo 



lm[(a xx -na xx )(a* xx - (77* + C)a* xx )}, 
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